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We propose an error correction coding algorithm for continuous quantum variables. We use this algo- 
rithm to construct a highly efficient 5-wavepacket code which can correct arbitrary single wavepacket 
errors. We show that this class of continuous variable codes is robust against imprecision in the 
error syndromes. A potential implemetation of the scheme is presented. 
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Quantum computers hold the promise for efficiently 
factoring large integers Q. However, to do this beyond 
a most modest scale they will require quantum error cor- 
rection Q. The theory of quantum error correction is 
already well studied in two-level or spin-i systems (in 
terms of qubits or quantum bits) |^,||,||J|,g,0] . Some of 
these results have been generalized to higher-spin systems 
|p|,P| JlO|Jll[| . This work applies to discrete systems like the 
hyperfine levels in ions but is not suitable for systems 
with continuous spectra, such as unbound wavepack- 
ets. Simultaneously with this paper, Lloyd and Slotinc 
present the first treatment of a quantum error correction 
code for continuous quantum variables fL2| , demonstrat- 
ing a 9-wavepacket code in analogy with Shor's 9-qubit 
coding scheme || . Such codes hold exciting prospects for 
the complete manipulation of quantum systems, includ- 
ing both discrete and continuous degrees-of-freedom, in 
the presence of inevitable noise jl3| . 

In this letter we consider a highly efficient and com- 
pact error correction coding algorithm for continuous 
quantum variables. As an example, we construct a 
5-wavepacket code which can correct arbitrary single- 
wavepacket errors. We show that such continuous vari- 
able codes are robust against imprecision in the er- 
ror syndromes and discuss potential implementation of 
the scheme. This paper is restricted to 1-dimensional 
wavepackets which might represent the wave function of 
a non-relativistic 1-dimensional particle or the state of a 
single polarization of a transverse mode of electromag- 
netic radiation. We shall henceforth refer to such de- 
scriptions by the generic term wavepackets @ . 

Rather than starting from scratch we shall use some of 
the theory that has already been given for error correc- 
tion on qubits. In particular, Steane has noted that the 
Hadamard transform 
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maps phase-flips into bit-flips and can therefore be used 
to form a class of quantum error correction codes that 
consist of a pair of classical codes, one for each type of 
'flip' ||. This mapping between phase and amplitude 
bases is achieved with a rotation about the y-axis by 
7r/2 radians in the Bloch sphere representation of the 
state. In analogy, the position and momentum bases of a 
continuous quantum state may be transformed into each 



other by tt/2 rotations in phase-space. This transition is 
implemented by substituting the Hadamard rotation in 
the Bloch sphere by a Fourier transform between position 
and momentum in phase-space. This suggests that we 
could develop the analogous quantum error correction 
codes for continuous systems jn|. 

We shall find it convenient to use a units-free notation 
where 



position — xx (scale length) 
momentum — p / (scale length) , 



(2) 



where a; is a scaled length, p is a scaled momentum and 
we have taken h = 1. (We henceforth drop the modifier 
'scaled'.) The position basis eigenstates \x) are normal- 
ized according to (x'\x) — 5(x' — x) with the momentum 
basis given by 
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To avoid confusion we shall work in the position basis 
throughout and so define the Fourier transform as an 
active operation on a state by 



T\x\ 
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where both x and y are variables in the position basis. 
Note that Eqs. (||) and (||) correspond to a change of 
representation and a physical change of the state respec- 
tively. 

In addition to the Fourier transform we shall require 
an analog to the bit-wise exclusive-OR (XOR) gate for 
continuous variables. The XOR gate has many interpre- 
tations including controlled-NOT gate, addition modulo 
2 and parity associated with it. Of these interpretations 
the natural generalization to continuous variables is ad- 
dition without a cyclic condition. That is, we take 



\x,y) = \x,x + y) 



(5) 



By removing the cyclic structure of the XOR gate we have 
produced a gate which is no longer its own inverse. Thus, 
in addition to the Fourier transform and this generalized 
XOR gate we include their inverses to our list of useful 
gates. This generalized XOR operation performs transla- 
tions over the entire real line, which are related to the in- 
finite additive group on R. The characters x of this group 
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satisfy the multiplicative property x( x + y) — x{ x )x{y) 
for all x, y € K and obey the sum rule 



(6) 



where x( a; ) = & iix - Interestingly, this sum rule has the 
same form as that found by Chau in higher-spin codes 
JTo| . Once we have recognized the parallel, it is sufficient 
to take the code of a spin-i system as a basis for our 
continuous- variable code. 

Based on these parallel group properties, we are 
tempted to speculate a much more general and funda- 
mental relation: We conjecture that n-qubit error cor- 
rection codes can be paralleled with n-wavepacket codes 
by replacing the discrete- variable operations (Hadamard 
transform and XOR gate) by their continuous-variable 
analogs (Fourier transform, generalized-XOR and their 
inverses). As a last remark before embarking on the nec- 
essary substitutions (in a specific example), we point out 
that the substitution conjecture is only valid for qubit 
codes whose circuits involve only these (H and XOR) el- 
ements. We shall therefore restrict our attention to this 
class of codes. 

An example of a suitable 5-qubit code was given by 
Laflamme et al. |{L6| . We show an equivalent circuit in 
Fig. [j] |l7]]. As we perform the substitutions, we must 
determine which qubit- XOR gates to replace with the 
generalized-XOR and which with its inverse. To resolve 
this ambiguity, two conditions are imposed. First, we de- 
mand that the code retain its properties under the par- 
ity operation (on each wavepacket). We conclude that 
either gate may be chosen for the first operation on ini- 
tially zero-position eigenstates. Ambiguity remains for 
the last four XOR substitutions. As a second step, the 
necessary and sufficient condition for quantum error cor- 
rection [U||: 

^encode ^encode) = S ( x> ~ x) X a p , ^ a, , (7) 

must be met. Here l^encode) encodes a single 
wavepacket 's position eigenstate in a multi- wavepacket 
state, £ a is a possible error that can be handled by the 
code and \ a p is a complex constant independent of the 
encoded states. [Condition (H) says that correctable er- 
rors do not mask the orthogonality of encoded states.] 

In the case of a single wavepacket error, for our 5- 
wavepacket code, it turns out that amongst the condi- 
tions of Eq. (fl) only ( x ' ^L^\ x enco d e ), havin S 
errors on wavepackets 4 and 5, is affected by the ambigu- 
ity (see detail below). An explicit calculation of all the 
conditions shows that the circuit of Fig. || yields a satis- 
factory quantum error correction code (as do variations 
of this circuit due to the extra freedom with respect to 
the choice of operator acting on wavepackets 1-3). By 
analogy with the results for higher-spin codes, we know 
that this code is optimal (though not perfect) and that no 
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FIG. 1. Quantum error correction circuit from |l7t| . The 
qubit \ip) is rotated into a 5-particle subspace by the unitary 
operations represented by the operations shown in this circuit. 
Note that the 3-qubit gates are simply pairs of XORs. 



four- wavepacket code would suffice |L(J. The code thus 
constructed has the form 
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c encode/ 
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dwdydze 2i{wy+xz) 

x \z, y + x, w + x,w — z,y — z) . (8) 
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FIG. 2. This 'circuit' unitarily maps a 1-dimensional sin- 
gle- wavepacket state \ifj) into a 5- wavepacket error correc- 
tion code. Here the auxiliary wavepackets |0) are initially 
zero-position eigenstates. For degrees-of-freedom larger than 
qubits the ideal XOR is not its own inverse; here the daggers 
on the XOR gates represent the inverse operation. 

Let us demonstrate the calculation of one of the con- 
ditions specified by Eq. ((7|): 



^encode I ^4a ^5/3|a; enC0 He) 



(9) 



dw' dy' dz' dw dy dz e ^v+^ w 'v' -x' z') 

xS(z' — z) S(y' — y + x' — x) S(w' — w + x' — x) 
x(w' - z'\£l\w - z)(y' - z%\y - z) 
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x (w — x' + x — z\£^\w — z) (y — x' + x — z\£p\y — z) . 

Making the replacements w — > w + z and y — ► y + z in 
this last expression we obtain 



3 — 2i(x' — x) 2 



dwdydz e 2i ^ x ' '-*)(*>+v+z) 



x(w - x' + x\£ f a \w}(y - x' + x\£p\y) 
S(x' - x) 



(10) 



dwdy (w\£l\w) (y\£p\y) = 5(x' - x) \ a p . 

For the other cases we find by explicit calculation, for 
wavepackets j ^ k, that 

Kncode^ ^l^encode) = S ( x ' - x ) X ^ • (11) 
For j — k this constant is found to be 

(12) 
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A a /3 = — I dw (w\£l £p\w) 



where C is formally infinite. 

We shall argue that this infinity vanishes when the syn- 
drome is read with only finite precision, which is always 
going to be the real situation. However, this requires us 
to demonstrate that our codes are robust: that for a suf- 
ficiently good precision we may correct single- wavepackct 
errors to any specified accuracy. In order to understand 
how the error syndromes are measured, let us consider 
a simpler code, namely, the continuous version of Shor's 
original 9-qubit code: 



1 



l^encode/ 



r 3/2 



dw dydz e 2lx( - w +y +z ^ 
x\w,w,w,y,y,y,z,z,z) , (13) 



where parity alone removes all ambiguity. (This code has 
been independently obtained by Lloyd and Slotine |fL2f.) 
Since this 9-wavepacket code corrects position errors and 
momentum errors separately, it is sufficient to study the 
subcode 



c encode) — \x,x,x) , 



(14) 



designed to correct position errors on a single wavepacket. 
The most general position error (on a single wavepacket) 
is given by some function of the momentum of that sys- 
tem £ (p) and need not be unitary on the code subspace 
[Eq. (ra)]. The action of such an error on a wavepacket 
may be written in the position basis as 

£&)\ x ) = \Jdydpe 2 ^y- x H{ P )\y) = j dy £{y)\x - y), 

(15) 



where £{x) is the Fourier transform of £{p). Thus the 
most general position error looks like a convolution of the 
wavepacket 's ket with some unknown (though not com- 
pletely arbitrary) function. Suppose this error occurs on 
wavepacket one in the repetition code (|l4]). Further, let 
us use auxiliary wavepackets (so-called ancillae) and com- 
pute the syndrome as shown in Fig. |^, then the resulting 
state may be written: 



dy £(y)\x-y,x,x,-y,0,y) 



(16) 
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FIG. 3. Syndrome calculation and measurement: A state 
with a single-wavepacket position error (here on wavepacket 
1) enters and the differences of each pair of positions is com- 
puted. The syndrome {si,S2,S3} may now be directly mea- 
sured in the position basis. 

Everything up till now has been unitary and assumed 
ideal. Now measure the syndrome: Ideally it would be 
{— y, 0, y} collapsing the wavepacket for a specific y. Cor- 
recting the error is now easy, because we know the lo- 
cation, value and sign of the error. Shifting the first 
wavepacket by the amount y retrieves the correctly en- 
coded state \x, x,x). Note that this procedure uses only 
very simple wavepacket-gates: The comparison stage is 
done classically , in contrast to the scheme of Lloyd and 
Slotine, where the comparison is performed at the am- 
plitude level and involves significantly more complicated 
interactions Jl2| ]. 

It is now easy to see what imprecise measurements of 
the syndromes will do. Suppose each measured value of a 
syndrome s'j is distributed randomly about the true value 
Sj according to the distribution p me as 

(s'j-Sj). We find 

two conditions for error-correction to proceed smoothly. 
First, Prneas^) mus t be narrow compared to any im- 
portant length scales in £(x). This guarantees that the 
chance for 'correcting' the wrong wavepacket is negligible 
and reduces the position-error operator to an uninterest- 
ing prefactor. If the original unencoded state had been 
J dx tp(x)\x) then after error correction we would obtain 
the mixed state 

dx' dxdz ip(x)ip*(x')p meas (z) 

x \x — z, x, x)(x' — z, x', x'\ . (17) 

Thus, unless Pmeas( x ) ^ s a ^ so narrow compared to any 
important length scales in ip(x), decoherence will appear 
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in the off-diagonal terms for wavepacket 1 of the cor- 
rected state ([l7]). This second condition is also seen in 
the quantum teleportation of continuous variables due to 
inaccuracies caused by measurement These condi- 

tions roughly match those described by Lloyd and Slotine 
p2[ . We note that any syndrome imprecision will de- 
grade the encoded states, though this precision may be 
improved by repeated measurements of the syndromes. 
For our 5-wavepacket example (||), syndromes consist of 
sums of two or more wavepacket positions or momenta 
and are measured similarly. 

It should be noted that Chau's higher-spin code JTcj ] 
could have been immediately taken over into a quantum 
error correction code for continuous quantum variables 
in accordance with our substitution procedure. How- 
ever, we have produced an equivalent code with a more 
efficient circuit prescription: Whereas Chau gives a pro- 
cedure for constructing his higher-spin code using 9 gen- 
eralized XOR operations, the circuit in Fig. [2] requires 
only 7 such gates or their inverses. In fact, we could run 
this substitution backwards to obtain a cleaner 5-particle 
higher-spin code based on Eq. (||). 

In order to consider potential implementations of the 
above code let us restrict our attention to a situa- 
tion where the wavepackets are sitting in background 
harmonic-oscillator potentials. By the virial theorem the 
form of a wavepacket in such a potential is preserved 
up to a trivial rotation in phase-space with time. The 
two operations required may be implemented simply as 
follows: The rotation in phase-space, Eq. (Q), may be ob- 
tained by delaying the phase of one wavepacket relative 
to the others, and the XOR operation, Eq. (^), should be 
implemented via a quantum non-demolition (QND) cou- 
pling. There exists extensive experimental literature on 
these operations both for optical fields and for trapped 
ions 

The conjecture put forth in this letter leads to a sim- 
ple, 2-step design of error correction codes for contin- 
uous quantum variables. According to this conjecture, 
any qubit code, whose circuit operations include only 
a specific Hadamard transformation, its inverse and the 
ideal XOR, may be translated to a continuous quantum- 
variable code, by substituting these operators with their 
continuous analogs and then imposing two criteria - par- 
ity invariance and the error-correction condition - which 
remove any ambiguities in the choice of operators. We 
demonstrate the success of this coding procedure in two 
examples (one based on Shor's 9-qubit code (2), and a 
second based on a variation of the Laflamme et al. 5- 
qubit code Jl6|,[l7]]). The 5-wavepacket code presented 
here is the optimal continuous encoding of a single 1- 
dimensional wavepacket that protects against arbitrary 
single-wavepacket errors. We show that this code (and 
in fact the entire class of codes derived in this manner) 
are robust against imprecision in the error syndromes. 
The potential implementation of the proposed class of 



circuits in optical-field and ion-trap set-ups is an addi- 
tional incentive for further investigation of the robust 
manipulation of continuous quantum variables. 
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